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Abstract 

Cyclic dominance of three species is a commonly occurring interaction dynamics, often denoted 
the rock-paper-scissors (RPS) game. Such type of interactions is known to promote species coexis- 
tence. Here, we generalize recent results of Reichenbach et al. (e.g. Nature 448, 1046 (2007)) of a 
four-state variant of RPS. We show that spiral formation takes place only without a conservation 
law for the total density. Nevertheless, in general fast diffusion can destroy species coexistence. 
We also generalize the four-state model to slightly varying reaction rates. This is shown both 
analytically and numerically not to change pattern formation, or the effective wave length of the 
spirals, and therefore does not alter the qualitative properties of the cross-over to extinction. 

PACS numbers: 87.23. Cc, 87.18.Hf, 02.50.Ey, 82.20.-w 
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I. INTRODUCTION 



Pattern formation and stability of ecological multi-species systems have attained a lot of 
■Merest recently flflQflfl. These two top.es are usuaUy .inked together by the well-known 
fact that spatial inhomogeneity can stabilize species co-existence |3|, |6[ , and understanding 
both the issues themselves and their connection is of great importance. 

An important class of model systems here are the so-called rock-paper-scissorsgames. 
They describe the dynamics of three species that cyclically dominate each other 



Typically, these systems involve a conservation law for the total density This can be 



imposed by writing rate equations such that the densities of the three species always sum 
up to unity, or employing lattice-based simulations where a site is always in exactly one of 
the three states. In this case, the system has a reactive fixed point which can be unstable, 
marginally stable, or stable depending on other properties of the particular setting, such 



as discreteness of the populations or dimensionality. There is a large 



very similar models with different microscopic update rules 12 



13 



14 



Dody of work on 



151 ]. and recently 



more complicated six-species systems with a similar conservation law have been studied 
16H . Several rock-paper-scissors-like processes have been identified in ecology, both spatial 
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181 ] and non-spatial [19|, |20|, as well as in other contexts not related to population 



dynamics such as the public goods game 



2l|. 



Recently, Reichenbach and co-workers have studied aversion of the rock-paper-scissors 

221 . In lattice simulations, this 



game in which the conservation law has been removed 0. 
is achieved by building the model on four states: the three original cyclically dominating 
states and a fourth one that denotes empty space. If diffusion is added, it does not interfere 
with the global conservation laws or absence thereof, but serves to set a length scale different 
from that set by the lattice constant. It has been shown that in this case the reactive fixed 
point is always unstable, and with diffusion in two dimensions spiral patterns form, similar 
to those in the complex Ginzburg-Landau equation (CGLE) l], or vortices in ecological 
systems 23, 24 1 . The main conclusion in Refs. 0, Q] has been that as a function of the 
diffusion constant D there is a cross-over from a reactive state with all three populations 
present to an absorbing state in which only one of the populations survives. This is argued 
to take place when the value of the diffusion constant is such that the spirals outgrow the 
system size. 
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Here we analyze these models further by considering both the four-state case and a 
three-state version with a conservation of the total density. We show that removing the 
conservation law gives rise to spiral formation that does not occur if the total density is 
conserved. However, in spite of this, there is a mechanism involving a diffusion- induced 
length scale that leads to a cross-over to an absorbing state. Therefore, conclusions regarding 
population stability are qualitatively the same for both cases. Second, all previous theoretical 
studies have assumed that the microscopic processes are rate-symmetric. In other words, 
one can cyclically permute the three populations without any change whatsoever. However, 



experiments can both explicitly show formation of spiral patterns 



and still be built such 



that the detailed pair-wise reaction mechanisms are qualitatively different for different pairs 



of species 



18] . It is apparent that a small rate-asymmetry does not change these properties. 



We introduce rate-asymmetry into the four-state model [5], and argue analytically, that the 
resulting system is still essentially the CGLE but only after an unimportant (in the first 
order) change of variables. This result, also confirmed by direct simulations, shows that 
already rate-symmetric theory and rate-asymmetric experiments are comparable. 

This paper is structured as follows. In Section II, we define both the three- and the 
four-state models with and without diffusion, and discuss in more detail what is previously 
known about them. In Section III, we present our results for all cases considered. Finally, 
Section IV ends the paper with a summary and discussion. 



II. THE MODELS AND EARLIER WORK 



The three-state RPS model describes cyclic dominance of three states, A, B, and C, and 
is defined by the following reaction equations and corresponding rates 

AB — > AA with rate fi 

BC -> BB with rate n (1) 
CA -> CC with rate p. 

In the mean-field approximation, the rate equations describing the evolution of the system 
are 



d t a = 


fiab 


— fiac 


d t b = 


[ibc 


— fiab 


d t c = 


jjbca 


- \xcb , 
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where a, b, and c are densities of the states A, B, and C, respectively. The rate equations 

have a reactive fixed point at a — b — c — |, which is known to be marginally stable. 

Such behavior is usually considered highly unrealistic in biologically motivated systems, and 

therefore the MF approximation of this system is considered to be of no relevance. 

The same model with discrete densities (i.e. finite populations) with noise has been studied 

□ 

in the mean-field like, non-spatial version [10j. In this case, the marginal stability of the 
fixed point remains, and the behavior of the system becomes a random walk in the space 
of marginally stable circular trajectories. Since the population is finite, this always leads to 
extinction in the infinite-time limit. 

However, in two-dimensional spatially extended systems the three-state RPS game is 
known to be stable, at least in large enough systems. The simplest known theoretical 
approach reproducing the behavior is the four-site approximation by Szabo and co-workers 



In empirical systems, the three states above are often considered to be three cyclically 
dominating strains of bacteria. In such settings, the dominance of strain A over strain B 
leads microscopically to deaths of individuals of strain B. In these cases, strain A does not 
reproduce immediately but non-occupied space is created, and the filling-in via reproduction 
of any strain is a separate process. Therefore, the following set of reaction equations has 
been proposed 

AB -» EA with rate a 
BC -> EB with rate a 

(3) 

CA — > EC with rate a 

XE -> XX with rate fi , 
where X can refer to any state and E denotes empty space. In contrast to Eqs. in which 
the only parameter /i sets the time scale, these reaction equations contain two independent 



parameters. The corresponding rate equations are [22] 

d t a = a[/i(l — p) — crc] 

d t b = b[fi(l - p) - aa] (4) 
d t c = c[/i(l - p) - ab] , 

where p = a + b + c is the total density. These equations have a reactive fixed point 



which is linearly unstable for all /i and a 22]. 



To generalize both models to their spatially extended versions, let the populations live 
on regular square lattices, the reactions take place only in nearest-neighbour contact, and 
amend the reaction equations with an exchange reaction 

XY -> YX with rate D , (6) 

where X and Y can denote any state (including empty space in the four-state model), and 
D is the diffusion constant. 

Exploiting the instability of the fixed point of Eq. Q a spatially extended version of the 
system has been reeent.y approximative.-, mapped Q to the eomp.ex Ginzburg-Landau 
equation (CGLE) In the mapping, one shifts the reactive fixed point to the origin, 
expands around it to find a two-dimensional invariant manifold of the dynamics, expresses 
the dynamics on the manifold, and performs a normal-form transformation, i.e. finds the 
quadratic transformation that is an identity mapping to first order and the result of which 
has no quadratic terms. Upon adding the diffusion according to Eq. ([6]), one arrives at a 
variant of the CGLE. 

In accordance with the known behavior of the CGLE, it has been found that the spatial 
four-state model with diffusion leads to formation of spirals [5]. They have a characteristic 
wave length that scales as the square root of the diffusion constant. When the wave length 
is of the order of the system size, there is essentially one spiral in the system. With even 
larger diffusion constants, the system behaves essentially as a completely coupled one and 
an extinction due to similar reasons as in the noisy non-spatial case takes place flOf . In 
other words, there is an absorbing-state cross-over as a function of the diffusion constant D 
when the spiral wave length reaches the system size. 

In direct numerical simulations, all variants of the model are simulated on a regular 
square lattice of size L x L. For each microscopic time step, a process (selection or diffusion 
in the three-state case; selection, reproduction, or diffusion in the four-state case) is chosen 
randomly with probabilities proportional to the rates, a random lattice site and its random 
neighbour are chosen uniformly at random, and the reaction is executed if allowed by the 
rules. The time is increased by At = 1/ (rL 2 ) where r is the sum of the rates for each case. 
The procedure is repeated until the time t reaches a predefined value. 
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III. RESULTS 



A. Three-state model 

1. Mapping to the CGLE 

To answer the question if explicit handling of empty space, i.e. using the four-state model 
instead of the more traditional three-state one, is really necessary to see the spiral pattern 
formation and the consequent absorbing-state cross-over, we map the three-state model to 
an equation resembling the CGLE as closely as possible. To start, set the timescale in 
Eqs. (|2j) by choosing ji = 1. Introduce new variables xa = a — ^, xb — b — |, xc = c — |, 
and use the conservation of the total density a + 6 + c= lto express the dynamics in terms 
of xa and as follows 

1 2 

d t x A = -x A + ^x B + x\ + 2x a x b (7) 
2 1 

d t x B = ~^ x a ~ -^x B - 2x A x B - x\ ■ (8) 

Note that the treatment here differs from that of the four-state model, since there is no need 
to reduce the number of dynamical variables by finding an invariant manifold. From Eq. (j7]), 
the linearization of the system around the fixed point is 



x B 

The eigenvalues of the linearization matrix are 




«- 1 XA \ = I I :• i i is) 



A = ±i^, (10) 

which recovers the known fact that the fixed point is marginally stable. 

To transform Eqs. (0) to a normal form, apply first the change of variables 

2 1 

VA = v^ XA + 7f XB ^ 

Vb = x B (12) 
to arrive at the rate equations, now expressed in terms of da and ys 

d t yA = ^yB + ^(y A -y 2 B) (13) 

dtVB = —^-VA - V3y A y B ■ (14) 



The reason for this particular change of variables will become clear when comparing the 
normal form to be calculated below to the CGLE. Namely, this transformation ensures that 
the linear terms are written in a directly comparable way. 

Now, let us turn Eqs. (fT3l) into their normal form. Here, the task is to find such a 
quadratic transformation from the variables ua and ys to a new pair of variables, say za and 
zb such that the y's and the z's coincide to linear order and that the rate equations for z do 
not contain quadratic terms. In general, finding such a transformation involves finding the 
inverse of a particular 6x6 matrix, and particular corner cases of singular matrices could 
cause the transform not to exist. However, in the present example, this does not happen. 

By imposing the restrictions above, we find the necessary transformation to be the fol- 
lowing 

za = Va + VaVb (15) 
zb = Vb + ^y 2 A ~ ^yl (16) 

and the resulting normal form 

d t z A = ujz b + c 2 z 2 A z B + c 2 z 3 B (17) 
d t z B = -ujz A - c 2 z\ - c 2 z A z 2 B , (18) 

AT AT 

where u = and c 2 = — with fourth-order and higher terms in i/a and ys omitted. To 
connect this to the CGLE, consider the complex variable z = za + izb- Now, Eqs. (I17p with 
added diffusion can be cast as 

d t z = iuz + DAz + ic 2 \z\ 2 z . (19) 

Further, changing to a rotating coordinate system (replace z by ze tuJt ) removes the purely 
imaginary linear term, and the resulting equation reads 

d t z = DAz + ic 2 \z\ 2 z , (20) 

which is already close to the standard form of the CGLE 

d t z = z + (1 + ib)Az - (1 + ic)\z\ 2 z. (21) 

Comparing Eqs. ( J20l and ( I2TT) reveals a crucial difference: in the CGLE there is a non-zero 
linear term whereas in the normal-form complex PDE the three-state model maps to, there 
is none. This is a direct consequence of the stability properties of the original MF equations. 
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To further argue why omitting the linear term matters, consider a variant of the CGLE 
with a more general linear term vz 

d~ t l = vz + (1 + ib)Az- (1 + ic)\z\ 2 z. (22) 

Given v > 0, Eq. (j2"Tl) is recovered after the scalings z = \fvz 7 t = 1/ut, and x = \fvx 
Substitute a generic single-spiral solution l| 

z = a (r)e <(ut± ^ (r)) , (23) 

where a(r) and ip(r) are real functions and r and are the spatial coordinates in the 
cylindrical coordinate system, with the spiral core at the origin, to Eq. fl22|) . Consider the 
limit r — > oo assuming that a(r) tends to a constant (i.e. the solutions are bounded) and that 
ip{r) tends to qr, i.e. asymptotically the spiral looks like a plane wave with wave number q. 
One finds that for the single-spiral solution to exist, the equality 

a(r) ~^a = - q 2 (24) 

has to hold. Since a(r) is real, spiral solutions with a finite wave length exist only for positive 
v. In other words, there are no spiral solutions in the three-state RPS game. The same 
argument can be repeated for plane waves as well, leading to exactly the same condition. 

We have also verified the conclusion by direct numerical simulations of the three-state 
model. Example configurations are plotted in Fig. [T] in panels (d), (e), and (f). There are 
indeed no signs of spiral formation, but a growing diffusion constant does change the spatial 
correlations in the system. Below, we show that this leads to an absorbing-state cross-over 
as a function of the diffusion constant regardless of the absence of spiral pattern formation. 



2. Extinction cross- over 



Both the three-state and the four-state models have an emergent length scale in the 
diffusive regime. In the four-state model, this is the spiral wave length, and in the three- 
state case it is simply the correlation length. Since the four-state model has an extinction 
cross-over as a function of the diffusion constant because of the length scale, it is natural to 
suspect that this is the case in the three-state model as well. We have performed extensive 
numerical simulations of both models to study the extinction probability as a function of 
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FIG. 1: (Color online) Example runs of the system in a lattice of size Lx L = 200. (a), (b) and (c): 
the four-state model with different diffusion constants and fj, = a = 1. (d), (e), (f): the three-state 
model with different diffusion constants and fj, = 1. (g), (h), (i): the rate-asymmetric four-state 
model with different asymmetries r = 1 + e, D = jj, = a = 1.0. 

the diffusion constants. The results are shown in Fig. [2] For the four-state model we recover 
the results of Ref. jlj] and the results for the three-state model show a similar cross-over. In 
both cases, the diffusion constant at the cross-over scales as D c ~ L 2 where L is the linear 
dimension of the lattice. 

Next, we argue that the extinction cross-over is caused by the presence of a correlation 



length, a pp 



diffusion 



25 



ymg 



26, 



an argument previously used in the context of the contact process with 



271 ] . Imagine the system in its steady-state, with a small spatially localized 
perturbation caused by noise. Due to the exchange reactions, this perturbation diffuses with 
diffusion constant 2D, where the factor of two comes from the fact that at each exchange 
reaction two particles moved. If the typical lifetime of such a perturbation is t^, it will 



diffuse up to distance x = y/2Dta before decaying. After estimating td, the cross-over should 
take place when L 2 = 2Dtd so that the properly scaled diffusion constant at the cross-over 
becomes 

D 1 
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FIG. 2: The extinction probability as a function of the scaled diffusion constant D/L 2 in various 
cases. For both the three- and the four-state model there is a cross-over to an absorbing state with 
two of the three subpopulations extinct. The location of this cross-over scales as D c oc L 2 in both 
cases and in the four-state case the introduction of rate-asymmetry does not affect the location of 
the cross-over for small asymmetries. 

To estimate the perturbation lifetime, we have computed the autocorrelation function of the 
time series of the densities. Fig. [3] shows the autocorrelation for several diffusion constants 
D and a corresponding time series. The immediate observations are that the autocorrelation 
function is independent of D and that it decays exponentially so that a well-defined timescale 
exists. From the linear fits, we can extract the time scale and we have found td ~ 66 (see 
Sec. II for the definition of unit of time in the simulations), from which we get that at the 
cross-over we have ~ 7.6 ■ 1CT 3 . This estimate corresponds well to the location of the 
cross-over in Fig. [2j 



B. Four-state model 



The previous formulations of both the three- and the four-state models have assumed 
invariance of the model under cyclic permutations of the three non-empty states. In par- 
ticular, this is seen in the rate equations (jSJ) and (jl]) where the reproduction rate /i is the 
same for all species. However, this assumption is not necessarily fulfilled by any realistic 
empirical scenario 18]. To find out if this matters, we define a rate-asymmetric variant of 
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FIG. 3: (a) The envelope of the autocorrelation function of the density of state A in the three-state 
model with different diffusion constants D. The symbols show the numerical results and the solid 
black lines are fits to the form a(t) oc e~*/* d . For a wide range of diffusion constants we find td = 66 
independent of D. (b) The corresponding time series of the densities of the three states show 
oscillations with wildly fluctuating amplitudes, often associated with almost-unstable dynamics 
and time-scale separation [281 ] . The parameters are L = 200 and \i = 1.0. 

the four-state model by the following reaction equations 

AB -> AE with rate a 

BC -> BE with rate a 

CA — >■ CE with rate a 

AE -> AA with rate fj, (26) 

BE -> BB with rate p, 

CE — > CC with rate r/i 

XY -> FX with rate D, 

where r = 1 + e (e > 0) is assumed to be close to unity, i.e. e C 1, and X and Y can 
refer to any state (including the empty state E). This is the simplest possible extension of 
the original four-state model that incorporates rate-asymmetry. The corresponding MF rate 
equations are 

d t a = a[pi(l — p) — crc] 

d t b = %(1 - p) ~ era] (27) 
d t c = c[rfi(l — p) — crb] . 
These have the reactive fixed point 

a = c=- ^ ;b = - r 4 • (28) 

(2 + r)n + a' (2 + r)/i + a K J 
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To simplify further calculations, we would like to introduce such a set of three dynamical 
variables a, b, c so that the reactive fixed point satisfies a = b = c. To this end, make the 
transformation a = a, b = b/r, c = c to arrive at the rate equations 

d t a = a[p{l — p) — crc] 

8 t b = b[p(l - P) - era] (29) 
d t c = rc[p(l - p) - ab] , 

where p = a + rb + c, and the reactive fixed point is a = b = c = ^ 2 +r)^+a • ^ s ^ ne ^ rs ^ 
step of mapping these equations to partial differential equations comparable to the CGLE, 
linearize around the fixed point to arrive at the linearization matrix 



D 



^ p pr p + a ^ 

p + a pr p 

,2 



(30) 



y pr pr + or pr J 

where D = (2 + r)p + a. Now, we can make some quick observations of the consequences of 
r > 1. First, the effect of b via the reproduction terms is multiplied by r which is natural 
since at the steady-state the increased reproduction of C has to be balanced by an increased 
density of B which dominates over C and therefore regulates it. Second, all time derivatives 
of c are multiplied by r which can be understood so that the faster reproduction of C alters 
the intrinsic timescale. 

From here, one would go further by first calculating the effect of small e on the two- 
dimensional invariant manifold of the system by expressing the normal vector of the manifold 
at the fixed point as its value in the rate-symmetric case and a first order e-correction to 
it. From there, one can continue by expressing the dynamics on the e-corrected manifold 
with two dynamical variables, finding the corresponding normal form transformation and 
applying it, assuming small e whenever necessary. While this procedure is certainly possible 
by brute force (which we have certified by carrying it out), the resulting expressions tend 
to get rather heavy since eigensystems of 3 x 3-matrices are involved, for instance. They do 
not appear to be useful for gaining physical intuition, and thus we resort to a qualitative 
argument of what the result, i.e. the CGLE-like PDE, necessarily is. 



If r = 1 (e = 0), the resulting normal-form PDE is 22] 



d t z = c x z + DAz - c 2 (l + ic 3 )\z\ 2 z, (31) 
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where z = za + iz B and za and zb are the two dynamical variables of the system. Here, the 
linear and cubic terms result from the normal-form transformation, whereas the Laplacian 
term comes from adding diffusion on top of the MF treatment. Similarly, upon carrying out 
the procedure outlined above, the e-perturbed terms are the linear term and the cubic term. 
However, only the perturbations in the cubic term are relevant. This can be justified as 
follows. Given general linear terms with four e-dependent coefficients for a fixed e, one can 
apply a linear transformation that diagonalizes the linear part. If the resulting eigenvalues 
are not real, they can be made such by transforming to a rotating coordinate system by the 
transformation z i— > ze 1 ^, as we did above to arrive at Eq. ( f20l) . Then, either one of the 
variables can be scaled such that the diagonal terms are equal. By these tricks, the form 
of the linear term can always be cast to be as in Eq. (13T1) . The remaining coefficient may 
well be e-dependent, but this plays no role as to whether the resulting equation has the 
properties and symmetries of the CGLE. 

So, the most general leading order correction from rate-asymmetry on top of Eq. ( !3~T1) has 
to be of the form 

d t z = ... + e(diz\ + d 2 z\z B + d 3 ZAZ B + d A z B ) + ie(d 5 z A + d%z\z B + d 7 z A z B + d 8 z B ) , (32) 

where the ellipsis stands for the non-perturbed terms, and the d^s are coefficients that in 
the most general setting are functions of the parameters fi and a of the four-state model. 
Now, to study the effect of the correction on spiral formation, let us express the CGLE and 
its corrected counterpart in polar coordinates. Write the phase-space coordinates as 



z = Re 1 ' 



■9 



(33) 



and the position-space coordinates as (r,<p). The CGLE is now 



d t R = aR + D(AR - R(V9) 2 ) - c 2 R 3 



(34) 



Rd t 9 = D(2V9 ■ VR + RA6) - c 2 c 3 R 3 . 



(35) 



And upon adding the nonsymmetric perturbation, we arrive at 



d t R = c x R + D(AR - R(V9) 2 ) - c 2 R 3 + eR 3 f(8) 



(36) 



Rd t 9 = D(2Vd ■ VR + RAO) - c 2 c 3 R 3 + eR 3 



9(0) , 



(37) 
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where / and g are smooth 27r-periodic functions with the property that their Fourier series 
do not have a constant term and have only a finite number of higher-order terms. 

The single-spiral solution of the non-perturbed CGLE is given by Eq. (1231) . To arrive at 
a similar solution for the nonsymmetric case, consider first the corresponding non-spatial 
dynamical system, i.e. Eqs. (}3~Tj) and (1321) with the spatial derivatives omitted. In this 
setting, both the symmetric and non-symmetric cases produce limit cycles. In the symmetric 
case, they are circles, and in the nonsymmetric case with small e one can envision them as 
perturbed circles. Then, the solution that is correct up to first order in e can be found by 
the change of variables that transforms the perturbed circles back to regular circles. Such 
change of variables in the first order in e is of the following form 

R = R(l + eh(0)) (38) 

with the function h(6) to be determined. Substituting this to Eq. (136]) gives 

= -F g 2 [l + eh(6)] 2 + F„(l - F 2 ) + eF*f(6) 

= -F g 2 + F (l - Fq) + e(-2F q 2 h(6) + F*f(9)) + 0(e 2 ) , (39) 

where the on the right-hand side (RHS) the e-independent part vanishes if z(f,<ft,t) is a 
single-spiral solution of the CGLE, and the whole RHS vanishes up to first order in e if we 
choose 

h(0) = 0/(0) • (40) 

As a conclusion, the perturbed system has a single-spiral solution where the oscillation 
amplitude has an additional phase-dependent prefactor that essentially cancels out the per- 
turbation to the circular form of the limit cycle of the corresponding non-spatial system. 
This does not change the wave length of the spirals (averaged over 9) which, in turn, means 
that the value of the diffusion constant at which the spirals outgrow the system size does 
not change, so that one expects that a small asymmetry does not change the location of the 
extinction crossover. 

We have verified this prediction also numerically. First, panels (g), (h), and (i) of Fig. [T] 
show direct simulations of the four-state model with increasing values of the asymmetry 
parameter r. One sees that the visual appearance of the patterns is modified due to the 
phase-dependent prefactor h{6) but that the effective wave length appears to undergo no 
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changes. Furthermore, we have systematically varied the diffusion constant D for different 
values of r studying the extinction probability. The results are shown in Fig. ([2]). The evident 
conclusion is that for small e the cross-over remains untouched. Finally, we have visualized 
the perturbation to the limit cycles of the non-spatial dynamical system in Fig. 




Z A 



FIG. 4: (Color online) The phase-space flow and limit cycle trajectories for the non-spatial rate- 
asymmetric four-state model with r = 1.0 (on the left) and r = 1.06 (on the right). The parameters 
are ji = a = 1. 



IV. DISCUSSION 

In this work, we have addressed two questions related to variants of the rock-paper- 
scissors game in which the total density is not conserved. First, we have shown that the 
spiral formation observed in the four-state case is a property of the four-state case only, 
and does not generalize to the more widely studied three-state rock-paper-scissors game 



with conserved total density. The spiral formation has been previously argued [5| to lead 
to a cross-over from coexistence of all three populations to an absorbing state with one 
surviving population and the other two extinct. The mechanism has been considered to be 
explicitly spiral-formation induced since the cross-over takes place at exactly that value of 
the parameters where the spiral wave length becomes equal to the linear size of the system. 
Here, we have shown that a similar cross-over takes place in the three-state model as well 
even though there are no signs of spiral (or any other visible pattern) formation at all. We 
have identified that the mechanism of the said cross-over is the appearance of a length scale 
from the interplay between fluctuations and diffusion. 

This result has direct consequences regarding the stability of such systems or, put in 
other words, the extent to which biodiversity is maintained. Namely, the message is that 
in cyclically dominating systems of three species, fast diffusion or big mobility can destroy 
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species coexistence regardless of whether there is a conservation law of total density in the 
system or not. The mechanism can either show up as visible pattern formation (here, the 
four-state case) or not (the three-state case). However, attention has to be paid to the 
details of the system if quantitative predictions are to be made. The value of the diffusion 
constant at which the destruction of coexistence is observed can differ by more than an 
order of magnitude for different cases (see Fig. [2]). The result also hints to the direction 
that the law of conservation of total density might play a role in pattern formation in more 
complicated well (see Ref. [16), for instance). 

Furthermore, the previous studies have considered only the cases where the reaction rates 
between all pairs of species are equal. This could potentially be a serious limitation since 
such equalities do not exist in reality, generally [18] . We have extended the previously defined 
four-state model [5j to cases where the reaction rates are not equal. By looking at the first- 
order perturbation from the case with equal densities, we have been able to map this case to 
a partial differential equation that, in turn, can be turned into the complex Ginzburg-Landau 
equation with a transformation of variables that does not alter the qualitative properties of 
the system in the first order. In particular, the effective wave length of the spirals equals 
that of the case with equal densities. As a conclusion, a small asymmetry in the reaction 
rates does not change either the properties of the pattern formation nor the location of the 
cross-over to the absorbing state. Altogether, this result serves to give a partial explanation 
why pattern formation and cross-over to extinction are visible in experiments as well. The 
behavior at larger asymmetries remains to be explored. 
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